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Abstract. In this paper we analyze the motion of a network of three planar curves with 
a speed proportional to the curvature of the arcs, having perpendicular intersections with 
the outer boundary and a common intersection at a triple junction. As a main result we 
show that a linear stability criterion due to Ikota and Yanagida [13] is also sufficient for 
nonlinear stability. We also prove local and global existence of classical smooth solutions 
as well as various energy estimates. Finally, we prove exponential stabilization of an 
evolving network starting from the vicinity of a linearly stable stationary network. 

1 Introduction 

The motion of curves under the curvature flow has been widely studied in the past HH 
E] . Less is known about the evolution of networks under the curvature flow [U Q21 [20] • In 
this case the arcs in the network evolve in the normal direction with a speed proportional 
to the curvature of the arcs. At intersections with an outer boundary and at triple 
junctions boundary conditions have to hold. At the outer boundary one can prescribe the 
position (see [Ml HE]), or the angle with the outer boundary [HEE3]. At the triple junction 
Young's law, a force balance, leads to angle conditions. In this paper we are interested in 
the stability of stationary solutions to the curvature flow with a triple junction when we 
prescribe the natural angle condition of 90° at the outer boundary. For this case a linear 
stability criterion has been derived by Ikota and Yanagida [13] (see also p3]). We will 
demonstrate here that this criterion also leads to nonlinear stability. 

We now specify the problem in detail. Let tt be a bounded domain in M. 2 with C 3 - 
boundary dtt. We introduce a C 3 -function if> : M? — > K with Vip(x) ^ if ip(x) = such 
that 

tt = {x e R 2 | ij){x) < 0}, dtt = {x e R 2 | ip(x) = 0}. 

We search families of curves Tj, T 2 , and T 3 which are parameterized by time t and which 
are contained in tt. The three curves are supposed to meet at a triple junction p(t) 6 Q 
at their one end point and at the other end point they are required to intersect with dfl, 
see Figure [U We require for i = 1, 2, 3 

/3 i y i = 7 i « i on rj, (1.1) 

3 

^fr = at p(t), (1.2) 
r t ±dtt at Tin dtt. (1.3) 
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Figure 1: The curvature driven flow T t of a network with a triple junction at p(t) (left) 
and a steady state (right). 

Here V i and k % are the normal velocity and curvature of F\, respectively The constants 
(3 % and 7* are given physical parameters and T l are unit tangents to the curve which are 
chosen such that they point away from the triple junction. 

Equation (11.21) is a force balance and one can solve for the T*'s if the condition 

7* + 7 J < 7 fc for all {i,j, k} mutually different , 

is fulfilled. In the following we assume strict inequalities and an argument as in Bronsard 
and Reitich [I] gives that the angles 9 l between the tangents Tj and fulfill 

sin 8 1 sin 8 2 sin 6 3 

ryl ry2 ,-y3 

with < 9 % < 7r {i = 1, 2, 3) and 6 l + 6 2 + 8 3 = 2n. Existence of solutions to the evolution 
problem (ll.ll) - (ll.3p has been shown by Bronsard and Reitich [I]. We will show later that 
the energy functional 

3 

E[T t ] = J2l i m 

i=l 

where T t = U? =1 T l t and L\T)\ is the length of T], is a Ljapunov functional. The con- 
stants 7* can be interpreted as surface free energy densities (surface tensions) and the 
functional E is the total free energy of the systems. Sternberg and Zeimer [21] showed 
the existence of isolated local minimizers to E, which can be interpreted as solutions to a 
partitioning problem of two dimensional domains into three subdomains having (locally) 
least interfacial area. 

The paper is organized as follows. In the next section we present a way how to 
parameterize the problem. We derive a nonlinear nonlocal system of parabolic equations 
governing the evolution of curves driven by curvature. By means of the semi-group theory 
due to Lunardi [17J we prove local existence of a classical solution. Section [H] is devoted 
to the rigorous derivation of the linearized system of equations. We recall the result of 
Yanagida and Ikota stating an explicit condition for linearized stability of the governing 
system of equations. In Section @] we provide a usefull result guaranteeing local uniqueness 
of a stationary solution proved by the inverse function theorem and the result is to our 
knowledge the first result in this direction for networks. As a byproduct we also obtain 
an important bound for the displacement of the network in terms of the curvature. We 



Figure 2: Description of the local parametrization of the curve P. 



proceed by deriving useful geometric equations for the curvature and other geometric 
quantities in Section El Using the linearized stability criterion we show how to derive a 
priori estimates for Sobolev norms of the solution. These geometric equations are then 
used in order to prove usefull bounds for a solution. With the help of these energy type 
estimates we prove global existence of a classical solution. In the final section [7] we prove 
exponential stability of the stationary solution. 

2 Parameterization and local existence 

We consider line segments Tl, Tl and Tl meeting the outer boundary with an angle of 
90° at their one end point and having without loss of generality = (0, 0) T as their 
common other end point where we assume that ( 11.21) holds. Then we define an arc-length 
parameterization of (i = 1, 2, 3) as 

T: = ma)\ae[0,r]} 

with (0) = p* = (0, 0) T , G <9f2. In particular, we obtain that l l is the length of T\. 

Then we will extend as an arc-length parameterization of the full line which contains 
r*. We will now introduce a certain stretched coordinate system in order to allow for 
parameterizations of curves close to T\ (i = 1,2,3) over fixed intervals [0,P]. 

Let T\ be the unit tangent to Y\ pointing from the triple junction p* to the outer 
boundary and let iV* = RTl be a unit normal where R is the anticlockwise rotation by 
ti / 2. We then define 

// = (p,T;) R2 , riniq) = max{a | + qN\ G H} . 

We remark that the parameter fi l allows for a tangential movement of the triple junction 
along r*. We now set 

^( ( T,g, / i 4 ) = $;(e(a,g,^)) + giV:, 

where 

C(a, q, //) =n* + ^C4fi(?) - »') ■ 
Note that f (a, 0, 0) = a and «f (0, q, //) = /A 



We now define the parameterization of curves T = (T 1 , T 2 , T 3 ) close to r* = (TJ, r 2 , r 3 ) 
having their triple junction at the point p with the help of functions 

p l : [0, f] -> K 

which fulfill the conditions 

p l (0) = (p,iV:) R2 (f= 1,2,3) (2.1) 

(see Fig. ED- Set 

&(a) = V(a,p\a), l x l ), a E [0, f] . (2.2) 

Then the functions $ l parameterize the curves P in the neighborhood of as T* = 
{$ l ((x) | a G [0,Z*]}. Since $*(^*) = ^Tj, we have $*(()) = p*Tj +/o*(0)iV* = p, which 
implies that 

p 1 ^ 1 + p^OR 1 = p 2 T 2 + P 2 (0R 2 = p 3 T 3 + p 3 (0)iV 3 . (2.3) 
By virtue of the definition of p\ equation ( 12. ip , and Young's law 

3 3 

1=1 i=l 

we are led to 

3 3 

$>V = 0, £ 7 V(0) = 0. 

i=l i=l 

Furthermore, identities (12. 3p and the angle conditions give the following lemma. 
Lemma 2.1 Let us define the matrix 

. , c 3 c l s 2 s 3 c 3 ^ 1 

Q = L— C V C 1 C 2 S 3 s l 



1 — 



s 2 c 2 s 3 ctc^s 1 



where c l = cos9\ s z = sm9\ Then, for p = (p^p^p 3 ), and p = (p^p^p 3 ), it holds 
» T = Qp T (0). 

Proof. It follows from ( 12. 3p that 

p* = p J '(Ti,T:) R2 +fj(0)(Ni,Ti) M i 

for (i, j) = (1,2), (2,3), (3, 1). By the angle condition, we have 

(T;,Ti) R2 = cos£ fc , (7* N')** = cos(6 k + vr/2) = -sm6 k 

for {i,j, k) = (1, 2, 3), (2, 3, 1), (3, 1, 2). This implies p* - cV = -«V (0), so that we are 
led to 

-c 2 1 \ / p 1 \ / -s 2 \ / p : (0) 

















u 3 y 





1 - c 3 o p 2 = -s 3 p 2 (0) 

1 -c 1 / \ p 3 / \ -s 1 / V p 3 (0) 



Then we obtain 



/i 1 \ / cV 1 c 3 \ / s 2 \ / ^(0) \ / p\0) 

p 2 \=^-L-[ c 1 cV 1 ) ( s* ) ( p 2 (0) ] = Q| p 2 (0) 



p 3 / 1 1 c 2 c 2 c 3 / V s 1 / \ p 3 (0) J \ p 3 (0) 



which completes the proof. □ 
We now consider evolving curves 

V(t) = {^(a,t)\ae[0,r]}, 

where $*(■,*) are defined as in (TjOD such that p\-,t) (i = 1,2,3) satisfy (j2TT]) . We 
formulate the curvature flow for a network with the help of these parameterizations. For 
that purpose, the following quantities are needed 

Ti = T^lK, N i = T — 1 R$i T , 

I cr I I (J I 

r = ($j,iV% 2 , «*= ( J-T;,iV% 2 , z = l,2,3. (2.4) 
Then we obtain the following formulation for the curvature flow of a network: 

/3 J ($i,iV% 2 =Y(^-r;,iV\ 2 , 2 = 1,2,3, (2.5) 

3 

^ 7 y = at (7 = 0, (2.6) 

i=l 

(T\T 2 ) R2 = cos # 3 , (T 3 , T 1 )^ = cos^ 2 at <r = , (2.7) 
(iV\ V^($ 4 )) R 2 = at <r = f, * = 1,2,3. (2.8) 

Note that the conditions ip{¥{l\t)) = (i = 1,2,3) and ^(0,*) = $ 2 (0,t) = $ 3 (0,t) are 
always fulfilled with our choice of the parameterizations. We can formulate the problem in 
terms of (p^p^p 3 ) and obtain a system of three second order parabolic equations where 
each equation is defined on a different spatial interval [0,f]. We obtain one boundary 
condition at a — V for the z-th equation and the three equations are coupled through the 
three boundary conditions at the triple junction. 

Let us derive the form of the nonlinear system for p*(cr, t) (i = 1,2,3). Set X % = [0, 1 1 } 
and Ql otl = r x (t ,ti\. Equations (12.5j) give 

p\ = L i {jtJ„it)d{itJ J„ t it)+k i {JJ„it)A for (a,t) G Qi >T , (2.9) 
where L l (p\ p l a , p l ) and A l (p\ p l a , p l ) are 



with the notation J 1 = |$^.|, and the curvature k 1 = K l (p l , p l a , p l aa , p l ) is represented as 



1 



(¥ q ,mi) M2 pi a + {2(¥ aq , R¥ a ) M , + (% a ,R%) R 2}pl 



{J i (p\p^pW 

By virtue of Lemma [2.11 and (12.91) . we have 
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QT°( p\ )=Q\r°M(p, Pa ^)] T° L\u*W{v?) 
p\ j \ L 3 {u 3 )k 3 {u 3 ) 



where Q is the matrix as in Lemma I2TT1 T° is the trace operator to a = 0, i.e. T°f = f 
and M(p, p a , p) is the matrix 

M(p, p CT , p) = Id - diag(A 1 (n 1 ), A 2 (u 2 ), A 3 (u 3 ))Q 

with the notation u % = (p*,p^.,//). 

Remark 2.2 T/ie matrix M(p, p a , p) is invertible provided that 



o-=0' 



Y7 sup + sup 1^(01) <^o 

~i X t&[o,T] te[o,T] J 



(2.10) 



for some 5 > 0. Indeed, we have det M(p, p CT , p) = d{— 1 + (c 2 — A 1 (w 1 )s 2 )(c 3 — 
A 2 ^ 2 )^)^ 1 - A 3 ^ 3 )^)}, where d = -1/(1 - c l c 2 c 3 ). Then, A*(0) = (z = 1,2,3,) 
zmp/y ffra* det M(0, 0, 0) = d(-l + c 1 c 2 c 3 ) = 1^0. Since det M(p, p a , p) zs continuous 
with respect to p, p a , and p, we can conclude that for e < 1 there exists a <5 > such 
that det M(p, p CT , p) > 1 — e > provided that h2.l0\) holds. 

As a consequence, we are led to the following nonlinear nonlocal partial differential equa- 
tions for p l (cr, t) {i = 1, 2, 3): 



Pi = a 



( P \ pi, + a v, pi, //) 4V p, t°p«, »)r°pi ( 



3=1 



+f(p\ d aP \ T»p, T p a , p) for (a, t) e Q^ T 

where a l (p\ p l a , p l ) = 7 l //3 l { J l (p\ p l a , p 1 )} 2 and a!( (T°p,T°p a , p) is the (i, j)-component 
of the matrix 



oi(T°p, T°p CT , P) = Q T u M(p, p CT , p) diag(T u a 1 (n 1 ), TV(u 2 ), TV(u 3 )) . 



-i 



Furthermore, f 1 is a smooth function in K which is evaluated at lower order terms. Then, 
recalling the boundary conditions (12.6p - (l2.8p and Lemma \2. 11 we have the following non- 
linear system: 

3 

pi = a\ P \ P i, pv*, + av, pi, p i ) E a ^ T °P' T °p^ ri T °ri° 

3=1 

+f(p\ pi, T°p, T°p a , fi) for (a, t) E Q ,t, 
£tV = 0, ^V 2 ) = 0, g a (u a ) = at * = 0, (2 ' U) 

8=1 

6^(p i ,/i l )p^ + ^(p l ,/i i ) = ata = r (2 = 1,2,3), 
k V T = QT°p T for te (0,T] 

where u lj = (p 1 , pP , p\, pi, p}, p?) (j = 2,3) and 

S 12 (u 12 ) = * 2 ) R2 + * 2 ) R2 p 2 + ^ ^2 )R2p i + ^ ^2 )R2p i p 2 

g l3 (u 13 ) = *i) Ra + (* 3 , *J) Ra pi + (M/J, v^) R2 p 3 + (*s f v^;) R2 p3 p i 
-J 3 (p 3 ,p^,p 3 )J 1 (p 1 ,pi,p 1 )cos^ 2 , 

Now we are ready to state a local existence result. 

Theorem 2.3 (Local existence) Let a E (0,1) and let us assume that p E C 2+a (T l ) 
and Pq (i = 1,2,3) with sufficiently small norms HPoUc^^) an d \Po\ fulfill the compati- 
bility conditions 

3 

£ tVo = 0, £ 12 (<) = 0, 2 13 K 3 ) = at a = 0, 
b m(Po, Po)Ph,a + 9m(pl Ph) = a* = I* (* = 1, 2, 3), 
where u 3 = (p x Q , p> Q , pl a , fP Q (T , p x Q , p 3 Q ) (j = 2,3). Then there exists a 

3 

To = T (l/EllPollc^)) >°> 

i=i 

To being an increasing function of its argument and such that the problem ( EL IV) with 
{p l {-, 0), p J (0)) = (p , p ) (i = 1,2,3) has a unique solution 

{p\p\p\p\p\p 3 ) E C^\Q^ o ) x C 2 +^(QQ x C*^ 1 ^) x [C'[0,T ]] 3 
satisfying A2.10\) . 



In order to prove Theorem 12 .31 by using a contraction principle, we need some prepara- 
tions which consist of three steps: 1) the linearization of (12.111) around the initial data; 2) 
the verification of the complementary conditions for the linearized system; 3) the deriva- 
tion of suitable a priori estimate for solutions of the linearized system. 

Step 1. Let us derive the linearization of (12.111) around the initial data p l G C 2+a (I l ) 
and p l {i = 1,2,3). First we define differential operators as 

A = diag(a 1 (^),a 2 (^) ;a 3( lt 3))52 ; 

A 1 = diag^K), A 2 (^),A3( w 3 ))ai(r o po;T o p o)T o d 2^ 



and also define, for given functions (p\jf) G C 2+ot > l (Qi T ) x C X [Q,T\ (i = 1,2,3), 
F\a,t) = {a^-a'K)fc 

+ ^{A^K (T°p, T%, Ji) - A'KK(T p , 7%,, » )}T% a 

+f(p\pi,T^,T%,-p), 

where uj, = (pj^, p^, p = (Pq,Po>Po)> = Oo^o>/4))> = (pS&.T?), P = 
(p 1 , p 2 , p 3 ), and 71 = (/I 1 , /Z 2 , p 3 ). Then, setting i* 1 = (F 1 , F 2 , F 3 ), we have the linearization 
of the differential equation given as 

d tP T = A oP T + A lP T + F T (a, t). 

Let us derive the linearization of the boundary conditions. For p[, G C 2+a (I % ) and p l 
{i = 1, 2, 3), we define differential operators as 

{(7\7 2 ,7 3 ) for k = l, 

(6 21 (T°<)^,6 22 (T°Oa ff ,0)for k = 2, 
(6 31 (TX 3 R, 0, b^{T°u^)d a ) for k = 3. 
Here the components are represented as follows: 

6 2 V 2 ) = (*J,* 2 ) R2 + «,^) R2 p 2 - {(**,*> + MV}^cos£ 3 , 
6 2 V 2 ) = (**,*> + (plvfaft - {(* 2 ,v^) R2 + |^ 2 | 2 p 2 }^cos^ 3 , 



J 2 (u 2 

b^) = (VP 3 ,^) R2 + (^ 3 ,^) R2p l _ {(* 3 ,*3 )m2 + |^| 2 p 3 }^|cOS^ 2 . 

Also, we define differential operators as 

(6^(^,^)^,0,0) for k = l, 

(BlUr%, pi)d a ) i=l!2 , 3 = { (0,b 2 m (T l2 plp 2 )d a ,0) for k = 2, 

(0,0,6 3 n (r' 3 p 3 ,p 3 )9 CT ) for fc = 3, 



where T l% (i — 1, 2, 3) is the trace operator onto a = l l , i.e. T l *f = f\ t . Then we set 



B (0-M = {B ki (r p ,T p 0ja ,fj, )d a ) k ^ 1;2t3 , 

and also set, for {f,J?) E C 2+a ^(^) x C x [0,T] (i = 1,2,3), 
G 1 (t) = 0, 



g*(*) = r° v\u^) P l a + &»'Kv , CT - ^(^Kp 1 - pi) - ±-{uW - pi) 



dpi 



dg lj 



dp 



T ( U v)(pi-pi)-^(ui j )tf-pi) 



(D'g 1 ^^ + (1 - vW^iu 1 * - ttf ), - ttf > dt7 (j = 2, 3) 



Ah 



1,2,3), 



where Dg l i is the Frechet derivative of g lj and the bracket (-, ■) is the respective inner 
product. Then we have the linearization of the boundary conditions: 

B (0 ; d a )p T = G T (t), B (P; d a )p T = G T m (t) 

for G(t) = (G\t),G 2 (t),G 3 (t)) and G dQ (t) = (G^f), & m (t) , G 3 dU (t)) . 

Step 2. Let us verify that the complementary conditions hold for the linearized system. 
We refer to Lunardi [17] for more information on the role of the complementary conditions. 
For that purpose, we make some preparations. Let C (r, i() = (clq )ij=i 1 2,3 where 

rC 2 + r, a% =0 for i ^ j. 



Then we have 



det £ n 



n 



rC 2 + r 



Setting £ = (So) = (det £o) (£o) 1 ) we are led to 

3 



n 



rC 2 + r 



for i 7^ j. 



U^K)} 2 

The matrix of the boundary conditions at o = is denoted by 

(l\l 2 ,l 3 ) for fc = 1, 

(i6 21 (T°ni 2 )C, i6 22 (T°OC, 0) for fc = 2, 
(i6 31 (T°ni 3 )C, 0, i6 33 (T°ni 3 )C) for jfc = 3, 



W(0;i()) != i, 2 ,3 = 

and ah a = T by 

(i^(f;iC)) l= i, 2 ,; 



(i6jn(r'pJ,^)C,0,0) for k = 1, 
(0,i6g o (T^g, A ig)C,0) for fc = 2, 
(0,0,i6 3 Q (T i3 p 3 ,/i 3 )C) for fc = 3. 



To verify the complementary condition we will show that the rows of the matrix BqCq are 
linearly independent for all r £ C, Rer > 0, modulo the polynomial 



P(r,0 = i[{(-a(r)} 



where Co( r ) = ^(no)! 7 *! 56 + (i = 1, 2, 3). Here we note that Co( r ) are the roots of the 
polynomial det£o( r ;iC) which have a positive imaginary part. 

First let us verify the complementary condition at o = 0. To determine whether or 
not the complementary condition is satisfied, we have to verify that the system 

3 3 

]T^(0;iCR(C) = mod P(r,C) = n(C-a) (< = 1,2,3) 

k=l i=l 

has the unique solution (u\, u>2, oj^) t = or equivalently that U2, oj^) t = is the 
only vector satisfying Y^t=i u k&o % (0 ; i£) = mod C ~ Co (* = 1)2,3). That is, we 
may investigate that (ui, u>2, ^z) T = is the only vector satisfying Y^k=i u kBo l (® j iCo) = 
(i — 1, 2, 3). Thus it suffices to show that 



/ 



det 



7* 



7' 



7° 




^ 6 31 (T 



0„.13> 



^0. 



[2.12) 



o 



b 33 {T°ul 3 ) J 



Indeed, in the case (p , /i ) = (0, 0) (i = 1, 2, 3), we have 



det 



( 71 
6 21 (0) 

v b 31 (0) 



r 

6 22 (0) 




7° 




\ 



/ 



det 



6 33 (0) ) 



r 



r 




7 

sin 6 13 — sin 9 3 
-sin6> 2 sin6> 2 y 
-(7 1 + 7 2 + 7 3 ) sin # 2 sin 9 3 ^ 0. 



V 



Since the determinant in (12.121) is continuous with respect to p l Q and fj, l (i = 1,2,3), we 
are led to (12.121) provided that HPoHc 1 ^') an d Ia*qI are sma U enough. 

Next let us verify the complementary conditions at a — /A Similarly as in the case 
cr = 0, it suffices to show that 

diag(&k(T%^ * 0. 

Indeed, in the case (p , fi l Q ) = (0, 0) (i = 1, 2, 3), we have 

det{diag(6^(0,0),6L(0,0),6L(0,0))} = det {diag(|V^($I)|, |V^($ 2 )|, |V^(<f 3 )|)} 

3 

8=1 

Since the determinant is also continuous with respect to p l and p l (i = 1, 2, 3), we conclude 
det{diag(&^(T^/4),& 2 n (T^ ^ provided that ||p llc;i(2*) and 

are small enough. 



Step 3. Let us analyze the linearized system. Set X = C{X 1 ) x C(T 2 ) x C(X 3 ) and 
Y = C 2 (X 1 ) x C 2 (J 2 ) x C 2 (l 3 ). Define the realization of Aq in X with homogeneous 
boundary conditions as follows 

D(A ) ={<peY\<p,Ao<pe X, B (0 ; d a )<p = 0, Boil'; d^ip = 0}, 
A)¥> = A <p. 

Then we have the following lemma which, in particular, characterizes the interpolation 
spaces Da (P, oo). For a definition of Da (P, oo) we refer to Lunardi [Tj 



Lemma 2.4 

(i) The linear operator Aq : D(Aq) — > X z's sectorial. 

(ii) The characterization of the interpolation spaces Da {ol, oo) is given as 

^,oo) = /^ e ^]^^ = > . <f ^(0 1/2), 

\ {¥>eX i9 |Bo(0;9 ff ) y) = 0,Bo(/ t ;9a)v = 0} i//3e (1/2,1), 

w/iere X^ := C^T 1 ) x C 2/3 (J 2 ) x C 2/3 (J 3 ). 

Proof. See [221 Section 2] or adapt the argument in (TTJ, Section 3.1.5] for the case of 
systems with the estimates in j2j Theorem 12.2]. □ 

Set A\Lp = Ai^p for ip 6 D(Ao). Then we obtain the following lemma. 

Lemma 2.5 Lei A = Aq + Ai. Then A : D(A ) — > X is a sectorial operator. 

Proof. According to [T71 Proposition 2.4.1(h)], A is a sectorial operator in X if Ax is a 
bounded linear operator from D(A Q ) to Da (o/2, oo). Indeed, by means of (I2.13P and the 
definition of oi, we have, for <£> e D(Ao), 

P^II^^Aoo) < C ^||A J '(^)||c-(x^ HaxC^Po, ^°Po ;CT , M ) IU a/2 H^^iU 

< CoII^IId^o) 

where a G (0,1) and Co, Cq are constants which depend on UpoIIc^cp) an d | /j-q | . This 
completes the proof. □ 

Using an estimate as in the proof of [T71 Proposition 2.4.1(h)], we see that A : D(A) = 
D(Aq) — > X is a sectorial operator such that ci||^j||d(^ ) < < C2||^||r»(A ) f° r <P £ 

D(Ao) and some constants c±, c<i > 0. Hence Da(cx, oo) = Da (oc, oo) with equivalence of 
the respective norms. 

By virtue of Lemma [2.51 we find that A generates the analytic semigroup e tA . Then 
we are led to the following proposition guaranteeing the existence of a unique solution for 
our linearized system. 



Proposition 2.6 Let us assume that p G C 2+a (X l ) (i = 1,2,3) satisfy the compatibility 
conditions 

B (0 ; d a )p T Q = G T (0), B (f; d a )f% = G T m (0). 



For (^,7?) G ^ 2+a,1 (Q^ T ) x ^[O.T] (i = 1,2,3;, toe Unearned system 

d t p T = Ap T + F T (a,t), 

B (0 ■ d a )p T = G T (t), Boil'; d a )p T = G^(t), (2.14) 
p o (;0)=p o (z = 1,2,3) 

with the notation A = Ao + A± has a unique solution such that 

3 

i=i 

3 

2=1 

Proof. Adapt the argument in the proof of [TTI Theorem 5.1.19] to our linearized system 

dsn. □ 

Now we are ready to prove Theorem 12.31 by using the contraction principle. 
Proof of Theorem \2.3[ Set 

v = {( p ,m) g c 2 +^("g~) x c^HQ^) x C 2 +^(^) x [c^t]] 3 I 

(p l (-,0), / i l (0)) = (pj ) ,/ii ) ) (i= 1,2,3), WpW^a + Mci < K) 
for some bounded positive parameters K and T where 

3 3 

i=i i=i 
Then, for (p, p) G X>, we define the mapping 

T : X> 3 (p, p) i-> (p, p) 

where p is the solution of (I2.14p and p is given by p T = Q T°p T for such solution p. Once 
we prove that the mapping JF is a contraction on T> for suitable K and T, the mapping 
T has a unique fixed point in T> which implies that the nonlinear problem (12. lip admits 
a unique solution in [0, T]. 

Let us first prove that T maps T> into itself. Note that the lower order terms in F l 
and Gg n can be rewritten as 

f(u) = f(u ) + f (Df ( V u + (l- V )u Q ),u-u )d V , 
Jo 

9dn{<) = <?L«o) + / (Dgan(v< + (1 - vM,o)M ~ <o> d V- 

Jo 



Here u = (p\ ft, T°p, T°p ff , n),u = (p , p 0)(T , T°p , ^°Po,<x> Mo). «i = (pS M*), and u* 10 = 
(p ,p ). Moreover, Df 1 and -Dg^ are the Frechet derivative of /* and g l dn , respectively, 
and the bracket (•, •) is the respective inner product. Then, by means of (12.151) and 
p T = QT°p T , we have 

||p|| C 2 +a>1 + \\fx\\ c i < c{j2(\\Po\\c*+"(T*) + \\f(ua)\\c<*(i*) + \9dn(T lt ui)\) 

^ i=i 

+ ^(|^ 1 (r°^)||r% (T | + |^(r°^)||rv 0i(7 |)) + c^ 

j=2 > 

for v = min{a/2, (1 — a)/2}. Thus, choosing 

^ i=i 

+ ih(\V 1 (T ul>)\\T pl a \ + |^(T°^)||TV^|)}, (2.16) 

3=2 ) 

we conclude that there exists a time 7\ > such that 

||p|| C 2+cM + \\n\\ c i <K for T < T v (2.17) 

That is, T maps T> into itself. 

Let us prove that the mapping T is a contraction. For (pi,Mi)> (p2;M~2) e ^ with 
T < Ti, let 

(Pl,Ml) =^(Pl,Pl), (P 2 )P2) =-7 r (P2,P2) 

be the solutions associated with the linearized problem (I2.14p . Then, applying a similar 
argument to jU pp. 373-375] with /x T = QT°p T , we are led to 

llPi - P2\\c 2 T +a ' 1 + II Mi - P2IIC1 - Cjf^dlPi - pJc^"' 1 + II Mi ~ M 2 ||ci) 

for v = min{a/2, (1 — a)/2}. Thus, T is a contraction for T < T 2 , which satisfies 
CkT% < 1/2. Consequently, choosing T* = min{Ti,T 2 }, we find that T has a unique 
fixed point in T> for T < T*, so that the nonlinear problem (12.111) has a unique solution in 
[0,T] with (12.171) for T < T*. Further, checking the details of the estimate for the linear 
system, we obtain for t 6 [0, T] 

3 

£IIp*MII Ci+«(2Ti) < "^0 + C K T U , 

i=l 

where mo depends on ||Po Ilc7 1+a (x i ) an d |/x |- Then, there exists a time T3 > such 
that m + Ck_T u < Ivtiq for T < T 3 . Thus, choosing T = minjT,,,, T 3 }, we have 
J^^ =1 ||p*(-, t) Hc-i+c^i) < 2m for t G [0, T] with T < T Q . It is possible to guaran- 
tee 2m < So for sufficiently small HPoUc^^) an d |p l> where <5 is as in (12 . 1 Of) . By 
Lemma |2~TI |p J (£)| is estimated by X^Li llp l ('^)llc°(x 4 ), so that |p l (t)| can be smaller than 
5q if ||Po Ila 1 + a (x i ) an d |p | are small enough. This completes the proof of Theorem 12.31 
□ 



3 Linearization 

In order to linearize the nonlinear system ( 12.111) around the stationary solution T* = 
Ui=i r*, we need to establish the following properties of at (q, = (0, 0). 

Lemma 3.1 For the parameterizations = 1,2,3, in Section^ the following prop- 
erties hold: 

(i) ^,0,0) = and ¥(a,q,0) = ®{{a^ m {q) /P) + qNi{arf aQ {q) /?). 

(ii) *i(a, 0, 0) = Tl, %(a, 0, 0) = NI, and %{a, 0, 0) = (1 - a/l*)^. 

(in) K a (a,0,0) = (0,0f, K q (°A0) = (0,0) r , and 0,0) = (-l/f)^. 

(iuj ^(tr, 0, 0) = (0, Of and % alx (a, 0, 0) = (0, Of. 

Proof. By the definition of (i) is obvious. Let us prove (ii). Differentiating \l/ l (cr, 0, 0) = 
$*(cr) with respect to a, we readily derive ^(cr, 0, 0) = T*(<r). Applying a similar argu- 
ment to |7], we obtain {p*Qn(q)}'\q=o = 0. Thus (i) gives \E^(cr, 0,0) = Nl(a). Further, by 
the definition of £ l , we have 

£*(<t, 0, 0) = 1 - a/J* . 
It follows from the definition of and the Frenet-Serret formulae that 

* i >, q, aO = e>, ?, - ^irae^, 9, ^)) = o, ^rae^, <?, //>)■ 

Putting (q,[i l ) = (0,0), the third property of (ii) is derived. Finally, by using (ii), we 
have (iii)-(iv). □ 

By virtue of Lemma l3~Tl we are led to the linearization of (12. lip around the stationary 
solution = |J 4 3 =1 r*. 

Proposition 3.2 The linearization of $2.11\) is given by 

f ^Pl^Ma for <7G(0,r), 

3 

7 V = at a = , . . 

4-f 3.1 

p\ = pl = pl at (7 = 0, 

k p CT + /i*p l = at o = l\ 

for i — 1, 2, 3, where hi is the curvature of dQ at the point (Z 4 ) G fl <9f2. 

We remark that (13. ip corresponds to the linearized problem which was derived in a formal 
way by Ikota and Yanagida |13j . 

Proof of Proposition \3. 6 A Applying the same argument as in [71 Section 3] , using Lemma 
13.11 and k\ = (i = 1,2,3) we obtain from equations ( 12.91) and the boundary conditions 
at a — l % , the first and fourth equations in ( 13. ip . Thus we only derive the third equation 
of ( 13.11) . To simplify the notation, we set 

j\p\^) = J(p\p^,//), ^V^mV) = <Ptfjy.J„v\iJ) (j = 2,3). 



Then it is easy to obtain 



£=0 ' 



where d J 1 (0,0) is the Frechet derivative of J 1 at (0,0). Recalling the definition of g 1 - 7 
(j = 2, 3) and using Lemma I3~T| we have 

^(0,0,0,0)[p\^/iV'] = -l(TlTi)^ - i(Ti,TiW + {T^Nfapj, 



+(Nl,Ti) R ,pl - (-1// - cos^ 



for (j, jfe) = (2, 3) or (3, 2), where dg lj (0 } 0, 0, 0) is the Frechet derivative of g lj at (0, 0, 0, 0). 
Since the angle conditions at a = give 

(T^TiV = cos# 3 , (T^iVfV = -sin# 3 , (i\^T 2 ) R2 = sintf 3 , 
(Ti,T 3 ) M2 = cos# 2 , (1?, JV?) Ra = sin# 2 , (iV^T 3 )^ = -sin# 2 , 

it follows that 

«^ 12 (0,0,0,0)[pV,/iV] = (-p'+p^sintf 3 , 
«^ 13 (0, 0, 0, 0) [p 1 , p 3 , p 1 , p 3 ] = (p 3 - pi) sin fl 2 

at (j = 0. Hence the linearization of ^ (p 1 , p 7 , p 1 , p- 7 ) = (j = 2,3) is 

(-p 2 + p^) sin 6 3 = 0, (p 3 - pi) sin # 2 = at a = 0, 

so that, by virtue of Qi G (0, tt) (j = 2, 3), we have 

Pa = Pi = Pi at o- = 0. 

This completes the proof. □ 

In [13] Ikota and Yanagida investigated linearized stability for the curvature flow with 
a triple junction ( 13. lj) . They derived a criterion according to which one can determine 
whether the stationary solution is linearly stable or unstable. In what follows, we recall 
their linearized stability criterion. The main result of [13, Theorem 1.1] is concerned with 
the analysis of the self-adjoint eigenvalue problem associated to the linearized system of 
equations (I3.ip . We now recall their linearized stability result. 

Theorem 3.3 The maximal eigenvalue of the linearized problem $3. i.e. f3 l = 0, 
i = 1,2,3, is negative and the stationary solution is linearly stable if and only if one of 
the following conditions is fulfilled: 

(a) either all hi, hi, and h\ are positive, 

(b ) or, at most one of them is non-positive, and they satisfy 

V (1 + fhDKhl + 7 2 (1 + ?hl)h\hl + 7 3 (1 + l 3 K)hlhl > . 



We will also need a variational characterization of the linearized stability property. To 
this end, let us introduce the bilinear form 

^^] = ^l l [j\^) 2 ds + K{^)X=i^ (3.2) 
for all ip G £(r*), where 

3 

S(T.) = {(<pW,<P*) G H'iOJ 1 ) x H\Q,f) x JT x (0, 1 3 ) | ^7V(0) = 0}. 

i=i 

This bilinear form was also considered in [13J. The following lemma is a simple conse- 
quence of the variational characterizations of the largest eigenvalue. 

Lemma 3.4 Let A be the maximal eigenvalue of the time independent linearized system 
/ TO) , i.e. we set ft = 0, i = 1,2,3. Then 

3 

i4^]>(-A)X>il^ll^ ( rt) for <pee(r.). 

1=1 

Remark 3.5 In order to simplify the presentation, we will henceforth consider only the 
case ft = 7*, for i = 1,2,3. It is worth to note that the linearized stability criterion 
is invariant with respect to the positive constants ft > for i = 1,2,3. As it should 
be obvious from all the energy type estimates to follow, the full nonlinear stability of the 
stationary solution will not be affected by a different choice of positive mobility constants 
ft > 0/ori= 1,2,3. 



4 Uniqueness of the stationary solution 

In this section we prove the uniqueness of the stationary solution in a small if 2 -neighborhood. 
The inverse mapping theorem also gives a bound on the i7 2 -norm of p = (p l ,p 2 ,p 3 ) in 
terms of the L 2 -norm of the curvature n = (k 1 , k 2 , k 3 ). 
To this end, let us introduce the function space 

3 

M = {(p\p 2 ,p 3 ) G H'iOJ 1 ) x H\0,l 2 ) x H 2 (0,l 3 ) I 5>V(0) =0}. 

8=1 

Then p via parametrization (12. 2p defines a neighboring triple junction configuration 
such that the end points lie on dQ. 

Theorem 4.1 Let be positive. Then there exists a H 2 -neighborhood of p = in M., 
such that p = is the only solution of the problem 



k* = 0, <(^,ri) = vr/2, (4.1) 
<(r i (t),T i (t)) = cos# fc for k G {1, 2, 3} mutually different. (4.2) 



Proof. The idea of the proof is to use the local inverse mapping theorem for the cur- 
vature operator with appropriate boundary conditions. The positivity of i* will ensure 
invertibility of the linearization. 

Using the notation ( 12. 4ft of Section [2] we obtain 

k\ P ) = oV, pi, Mp(o)M, + f{ P \ p\, m(p(o))) 

with a l (p\ p l a , /i(p(0))) = l/{J l (p l ,p^,/x(p(0)))} 2 , a smooth function and a linear 
mapping p from R 3 to IR 3 . The boundary conditions in (14.11) and (14. 2 p can be written as 

SW),P^),MP(0))) = for t = 1,2,3, 
g%p(0),p a (0),p(p(0))) = for * = 4,5 

where g 1 , . . . , g 5 are smooth functions. We define Af = L 2 (0, 1 1 ) x L 2 (0, I 2 ) x L 2 (0, 1 3 ) and 
observe that solving (14.11) and (14.21) is equivalent to finding a zero of the mapping 

F :B £ (0)^^x^ 5 , F(p) = (K l (p),K 2 (p),K 3 (p),g 1 ,...,g 5 ) 

where B £ (0) is a ball of a radius e > around zero in the space M.. Since H 2 is embedded 
in C l the mapping is well defined. Arguing similarly as in Section [3] we obtain 

DF(0) : M -> A/" x M 5 

p ^ (pL ) pL ) pL ) pi + ^V 1 ,Pa + ^V ) Pa + ^ 3 p 3 ) pi-p^Pa-Pa)- 

Similarly as in [13], since J* is positive we can conclude that DF(0) is injective and hence 
the Fredholm alternative gives that DF(0) is invertible. Now the local inverse mapping 
theorem (see e.g. [21]) gives that there is a neighborhood around such that only p = 
solves g3D and (IP]) . □ 

It is worthwhile noting that the mapping F analyzed in the proof of the above theo- 
rem is in fact a local diffeomorphism. Therefore its inverse mapping is locally Lipschitz 
continuous. Hence we have the following corollary. 

Corollary 4.2 There exist constants C, 5i > such that 

\\p\\m < C\\k\\ L 2 

provided that ||k||^2 < 5% and p G Ai fulfills F(p) = ( 0,0,0,0,0). 

In other words, by means of the above theorem and its corollary, we obtained a bound 
on the if -norm of the solution p in terms of the L -norm of the curvature 
in the vicinity of the stationary solution p = provided that p fulfills the boundary 
conditions. This useful observation will be used several times throughout the rest of the 
paper. Although we could take the standard L 2 -norm in the corollary above we choose 
the suitable weigthed L 2 -norm defined in (16.21) in the corollary as this will simplify the 
further presentation. 



5 Governing equations for the curvature and other 
geometric quantities 

In order to show the global existence and the nonlinear stability of solutions for which 
the bilinear form of [7j is positive, we apply an energy method similar to the one used 
in [5] and p3]. For such a method it is important to derive evolution equations for the 
curvature. 

Let s be the arc-length parameter along the phase boundary T t and let X be a smooth 
map such that X(-,t) is an arc-length parameterization of T t with 

Tl = {X\s,t)\se[0,r\t)]} 

where r % is smooth such that L\T\] = r l (t) which is the length of T l t . Let N l (= N l (s,t)) 
be the unit normal vector of T\. It can be written as 

K ' \ smtu l (s,t) J 

Then we have 

f Nl = -K ir P, Tl = <jN\ 
\ W t = -u\T\ Tl = ulN 1 

where T l is the unit tangent vector of T\ and k % is the curvature of T\. In addition, we 
define 

V i = (Xi,N i )w 1 v l = (Xt,T% 2 

and hence 

X\ = V i N i + v l T\ (5.2) 
Differentiating 05.21) with respect to s and using (I5.ip . we have 

XI = V'N 1 + V i N t 8 + v\T l + v { T l s (5.3) 
= (V; + kV)JV* + (-k¥ + vi)7* . (5.4) 

Lemma 5.1 Let X 1 be a smooth arc-length parameterization as above. Then 

J t = V l s + reV, v\ = k 1 V { . 

Proof. Since Xj s = and X l s = T\ it follows from (EI]) and (E2D that 

wiN* = (y; + «V)JV* + (-k'V* + v\)T\ 

Thus we obtain the desired results. □ 

As a consequence of Lemma 15.11 we have the following lemma. 

Lemma 5.2 Let X 1 be a smooth arc-length parameterization of T\ as above. Then the 
curvature k 1 (i = 1,2,3) satisfies the evolution equation: 



Proof. By uj 1 s = k 1 and Lemma 15. 1\ we obtain 

4 = u\ t = u\ s = {v: + K v) s = v: s + «x + k*v* = v: s + (*«) v + 

This completes the proof. □ 

By the assumption that T\ meets dfl at the one end point with the angle 7r/2, we have 

^{X\r\t),t)) = 0, (W>(X*),iV% 2 = at s = r\t). 

Differentiating the identity ip(X l (r l (t),t)) = with respect to t and taking into account 
the transversality condition (S7ip{X l ), N' 1 )^. = and the governing equation X\ = V l N l + 
v % T % we can derive the following lemma. 

Lemma 5.3 At the "point where Y\ (i = 1,2,3J meets the outer boundary dVL we have 

^ + (7-7 = for s = r\t). 

Next we derive corresponding boundary conditions at the triple junction point pit). 
It is assumed that phase boundaries T 1 , T 2 , and T 3 meet at the triple junction. Let 
i,j, k G {1,2,3} be mutually different. Let T k and 7 fc be the interface and the surface 
energy density between phases i, j. Following the arguments in Bronsard and Reitich [lj 
the angles 6 l (i = 1, 2, 3) of the phases at the triple junction point p{t) fulfill Young's law 

sin^ 1 sin^ 2 sin# 3 

— r~ = — ~ = — 3~ <- 5 - 5 ) 

✓y-L ry4 ryO 

(see [23J). Young's law can be expressed as a force balance in the following form 

3 3 

J2YT* = Y,l l N l = 0. (5.6) 

i=l i=l 

Let p(f) G M 2 denote a triple junction. At the triple junction the following boundary 
conditions hold: 

X\0,t) = X 2 (0,t) = X 3 (0,t) (= p{t)), (5.7) 
iXl,Xi) R2 = cos9 k at pit) (5.8) 

for (i, j, k) = (1, 2, 3), (2, 3, 1), (3, 1, 2). Then we obtain the following lemma. 

Lemma 5.4 At the triple junction pit) we have the following equality: 

3 



^ 7 V = 0. 



i=l 



Proof. Differentiating (15.71) with respect to t, we obtain 



For % = 1,2, 3, it holds 



d {t= X '= X " = X '- ^ 



This fact and Young's law imply that 

i=l i=l i=l 

Hence the proof is complete. □ 

In the next lemma we derive evolution equations and boundary conditions for the 
curvature. 

Lemma 5.5 A smooth solution of the curvature flow equations 

V i = k\ i = 1,2,3, (5.10) 

with the boundary conditions 

<(P(t),P'(t)) = cos0 fc for i,j,k e {1,2,3} mutually different, 
<(dn,Ti) = n/2 at V\ n dtt, (5.11) 

dv\ c on 

fulfills when expressed in the above arc-length parameterization the evolution equations 

4 = <s + (^) 3 + <v\ 2 = 1,2,3. 
Furthermore, at the triple junction p(t) we have 

3 

5>V = 0, (5.12) 

2=1 

K \ + kV = k 2 s + kV = k\ + «V (5.13) 

and at T l t fl dQ the identity 



K l s + h l K l = (5.14) 
holds. Here h 1 is the curvature of dQ at the point X l (r l (t),t) £ fl dQ. 

Proof. From f)5.9p we deduce 

(%N% 2 = (XIN%. = V\ 
This fact and Young's law imply that 

i=X i=l i=l 



Since V 1 = k 1 , we are led to (15.1 2\ . 

Differentiating (15.81) with respect to t, we obtain 



(Xlt, X J S ) R 2 + (XI, X J st ) R 2 = 0. 



It follows from (EU and (l5~TTD that 

(V; + «V) sin 6 k + (V> + reV)(- sin # fe ) = 0. 
By < 6 k < 7r we have sin # fe > 0, so that 

V l s + kV = V s j + K j V j . 
Finally, Lemma 15.21 and the fact that V 1 = k 1 imply (I5.13P . □ 

6 A priori estimates and global existence of a smooth 
solution 

The purpose of this section is to derive a priori estimates guaranteeing global existence of 
a smooth solution and its convergence to a steady state. First, we derive a priori estimates 
for the L 2 -norm of the curvature. Next we proceed with higher order energy estimates 
yielding a priori estimates for the if 2 -norm of the curvature. As a consequence of these 
estimates we will be able to prove exponential decay of the if 2 -norm of the curvature. We 
remark that due to the parabolic regularization property the solution of Theorem 12.31 will 
become smoother for positive time such that all derivatives in the following computations 
exist. 

6.1 First order a priori estimates 

Let us define the energy functional 



3 



i=l 



where T = U? =1 T l and L[T l ] is the length of T\ Further, throughout Section El we use 
the following notations: 



V^Wlp — / lv 9 T c ^ s ] (1 — V < °°)> Hv^lU 00 = esssup |<^ l |, (6.1) 





IMU< 




(6.2) 



for a vector function ip = ((p 1 , if 2 , ip 3 ). 



Lemma 6.1 A smooth solution of $5.10\) - [5.11\) fulfills the following energy type identities 



(i) jE[Vt] + Ml* = 

3 



where h % is evaluated at X*(r*(t),t). 



s=0 



Proof. By means of the identity L[rj] = r l (t), Lemmata 15. 1[ [5731 and l5~4"l we have 

3 3 

dt 



im) = EtV)' = -EtW,*w(o,*)} 



i=l i=l 

3 /t 3 /» 3 /» 

-EW «> = -E'y < / « i v*cfa = -E7 < / c**) 2 * 



1=1 1=1 i=l 

In order to prove (ii), we compute 



2 / ds = 2 + (k*) 3 + «y } tfa. (6.3) 



The identity 

d 



dt J r 

and Lemma [5.31 imply 



-77 / (n i ) 2 ds = 2 I ^Klds+i^fl^Jr^' 



I s=r' 



2/ K*Klds=j[ i (K i ) 2 ds+(K i ) 2 V i \ s=ri . (6.4) 



r; ut </r 



For the right hand side of equation (16.31) we can use the boundary condition (15.141) . 
Integration by parts yields 



= 2{[« < «J^- / (^ s ) 2 \ds + 2 [ (K i fds + 2 [ K^yds 

= -2\f (K i s ) 2 ds + h i (K*) 2 \ s= X-2K^ s \ s=0 + 2 [ {K i fds + 2f ^xiyds. (6.5) 
Using the identities v l s = k % V % = (ft*) 2 and integrating by parts we obtain 

Thus we have 



It follows from fH}, O, and USD that 



1 / H-fds = -2{ I (H'j'ds ■ h'(,')\ r \ ■ I (KTr/.s 



-2k'k*| „-(«*) VI .. 

s I s=0 v ' I s=0 

By (15121) and (15131) we have 

3 3 

ls=0 



E^{ 2 ^L=o + (^) v L=o} = 2 E7V(<+«v)| s=0 -^f(/, 

i=l i=l i=l 

3 



s=0' 



i=l 

Thus we are led to the identity 

3 3 



»=i i=i J i=i ^ r t 

3 



1 » ' ls=0 
i=l 

Since V 1 = k % the proof of the lemma follows. □ 
Let us define a bilinear form / as 



i\2| 



I s=r' 



i=1 Uo 
for y> G £(T) where 

3 

£{T) = {{<pW,<f?) e ^(O.r 1 ) x ^(0,r 2 ) x ^(0,^) | $>¥(0) = 0}. 

i=i 

Since V G £(T) we can rewrite the statement ii) of Lemma 16. II as 

±\\k\\1 + 2I[V,V\ = rt4 + ^7V)V| s=0 . (6.7) 

i=i 

The following lemmata are crucial in the derivation of a priori estimates. 



6.2 Higher order estimates for the curvature. 

We define the averaged curvature along the curve T l t as K l av = -^py J Tl n l ds. 

Lemma 6.2 The following estimates for a C 3 -curve T] hold true. 



V \\^-<v\\l- < (L\r$M\\i?, W av \ < -^—t\W\\l>, 



(Hi) There are m 1 = m,i(L[r\]) and m 2 = m 2 (l / L[Tl\) such that 

Proof. The estimates in (i) are established in a standard way and we do not present 
details here. By the estimates in (i), we have 



l«*IU~ < Ik - «iJU« + Kv\ < (£[ri]) 3 ||/4lU* 



It implies that 



nTds 



< ll^ll 2 Wk'W 2 

— II II L°° II II L 2 



< \(l\t>])*\\k*\\v + 



-\\ K IU 2 c \\ K Hl 2 



The statement (ii) now follows from the elementary inequality (a + b) 2 < 2(a 2 + b 2 ). 
Let r be such that K % (r l ,t) — K % av (i) = 0. We then obtain 



i K% K av) 3 \s-- 



=o = -3 / 
Jo 



K av) ds. 



This implies that 



|(/^-<J 3 Uo| < 3 / \Kl\\^-^ v \ 2 ds 



< 3H/C* - <JU-||4IU 2 II^ - <JU 2 

< 3(L[ri])i||<||| 2 {||^|| L2 + (L[ri])^|<j} 

< 6(L[rj])l||^|| L2 ||<||| 2 . 

Note that (V) 3 = (k* - <J 3 + 3^<„(^ - <J + (<J 3 . Then it follows that 

< L[rj]|<j||<||| 2 + (L[rj])§|<j 2 ||<|| L2 



K 'II2 2 + ll^slli 2 



= ^(i[r;])* + i}|| K -|u s ||<||l a + ^|| f c-||!, 



Thus we have 

K*i s =o) 3 | < \(K i -< v )%= \+nK i < v (K i -< v )h~ + K vi 



v av | 



2 j n ii* n ,ii^ [2L[rj] (L[rj]) 3 / 2 

Hence the proof is complete. □ 

Lemma 6.3 For smooth solutions of $5.10\) . ^5.11\) we have 

( V \ v\ vY = q(v\ v\ vy, 

for any t G (0, T], where the matrix Q was defined in Lemma \2. 1\ 
Proof. At the triple junction p(t), we have, for all i,j G {1, 2, 3}, 

^ = tfT + V i N i = v j T j + V j N j . 
dt 

Taking the inner product with T l we obtain 

v { =v j {T\T j ) + V j {T\N j ). 
By (EHD, (T\T j ) = cos6 k and (T\N j ) = -sin6 k . Thus we derive 

j _ c k v j = -s k V j . 

If we solve this with respect to (v 1 ,v 2 ,v 3 ) T , we are led to the desired result. □ 
By Lemma [6.31 and V 1 = k\ we have 



L ' 2 ' 



3 

7V) 3 Ln- (6-8) 



\s=0 



i=l M=l J M=l ' i=l 

6.3 Structural stability of the bilinear form 

The aim of this subsection is to show that positivity of the bilinear form 7* is invariant 
with respect to small perturbations of the curve parametrization p % . More precisely, we 
will show that the bilinear form / is positive definite provided that is positive definite 
and p = (p x ,p 2 ,p 3 ) is small in the C 1 -norm. Taking into account Corollary 14.21 and 
the continuity of embedding H 2 <^-> C 1 the positive definiteness of the bilinear form is 
preserved if the L 2 - norm of the curvature k, = (/c 1 , k 2 , k 3 ) is small. 

According to Lemma 13.44 7* is positive provided the maximal eigenvalue for the 
linearized problem is negative. The following lemma is a direct consequence of [TBI 
Lemma 3.1 and Prop. 3.3]. 

Lemma 6.4 

(i) Let A be the maximal eigenvalue of the linearized problem. For e > there exists a 
5 > such that, for any perturbation \h l — h l J\ < 5 and |L[r i ] — L[r*]| < 5, i — 1, 2, 3, 
we have 

I[ (p ,< fi } > (-\-e)\\<p\\ 2 L2 for<pe£(T). 



(ii) There exists a c > such that 

c \\<P.\\h < J [<?> y] + llvllia /° r V e £( r )- 



Using Lemma 16. 4[ we obtain the existence of constants 5* > and c* > such that 



Jfo¥>] > -yllvllia + c»ll¥».llia for ip G £{V t 
provided that, for i = 1, 2, 3, 

\ti-hi\ <5*, \L[r]-L[ri]\<6.. 



(6.9) 
(6.10) 



Lemma 6.5 We have the following estimates. 

(i) There exist constants 5 2 , C > such that \h l — h\\ < C||p l ||co(:r<) provided that 
1 1 Pile* (J 1 ) < $2- 

(ii) There exist constants 63, C > such that 

\L[r}-L[Tl]\ <C\\p\\ c i and \p\ < C||p|| c o 
provided that \\p\lc 1 < ^3? where p is the triple junction ofT = |J^ =1 T\ 
Here we have denoted ||p||c*+a = Y^h=i \\p l \\c k + a (i i ) f or & G N U {0}, a G [0, 1). 
Proof. To prove (i), we recall that kqq(X 1 ) is represented by 

1 



■([D^xyTaniX^TgniX*))*. 



(6.11) 



Since the right hand side does not depend on derivatives of p\ the mean value theorem 
implies the second inequality of (i). 

In order to prove (ii), we have to analyze properties of the function fi dn introduced in 
Section El From the definition it follows that fM dQ (0) = l l and p* + p, l an (q)Tl + qNl G d£l. 
Therefore ip(p* + ii dn (q)Tl + qN'l) = and so we can compute the derivative of p, Qn (q) as 



d 



/4n (?) 



Since (Vip, iV*) K 2 = on dfl we obtain J^/^q(0) = 0. Now, by taking the second derivative 
of 

^an an d taking into account the expression for the curvature hi at the intersection of 
dil and T'^ we obtain ^/x a n(0) = h\. Thus 

/4n(<z) = + ^*<? 2 + o(q 2 ) as g -> . 
We recall that the parameterization $ l of the curve T l is given by 

=p*+ [> + ^{/4 n (pV)) - A*'}] ^ + P>)^ • 

Using the above property of the function n an and the fact p T = Qp(0) T we obtain 



|L[r] - L[rj 



I) da 



<C\\p\ 



c 1 



Similarly, as p = p}Tl + //(0)iV* we obtain \p\ < C\\p\\co. With this all statements of the 
lemma have been shown. □ 



6.4 Exponential stabilization of the solution 

Lemma 6.6 Let A be the maximal eigenvalue of the linearized problem. Assume that X 
is negative. Then there exists a S3 > such that 

||^)|li 2 <e A '/ 2 || K (0)||| 2 for te[0,T], 

f \\^{r)\\l 2 dr<U^mh 
Jo c * 

provided that ||K(t)||| 2 < S4 holds on [0, T] where c* > is a constant as in ( Iff.gj) . 

Proof. According to Corollary 14. 2} there exists a C > such that ||p||ir a < C|I k IIl 2 f° r 
||k||^ 2 < Si. By ( 16. 7p we have 

j t \\K\\l 2 + 2I[V, V] = + $>V)V| s=0 . 

i=i 

Let us first choose S2 € (0,5*) H (0, Then, it follows from Lemmata 16.21 and 16.51 and 
the inequalities (16. 8p . (16. 9p . and (16.101) that there are C > such that 

j t \\ K \\l 2 + (-A)||V||f 2 + 2c.HV.Hi, < £7(11*11* + + KlliO- 

Since V 1 = k 1 , we are led to 

^ll^ll* + {(-a) - c(\\k\\v + HI*)}NI* 

+ |2c* - C7(||k|| l2 + ||K||ia) }||K.||ia < 0. (6.12) 
Then, we choose a constant 5 4 > satisfying 

< ft < minjl, g) ,^,(|j) 2 }. 
If we assume ||K(t)|| 2 2 < &i for t e [0, T] we then have 

+ + < 0. (6.13) 

Using the Gronwall inequality we obtain the desired result. □ 



6.5 Higher order energy inequalities 

So far we have shown the exponential decay of the L 2 -norm of the curvature «\ In order 
to prove stabilization of the curvature in the stronger C 1+Q -norm we need to derive higher 
order energy type inequalities. These estimates will enable us to conclude convergence 
of the curvature to zero in the C 1+Q -norm. In order to derive higher order estimates 
we differentiate the curvature equation (see Lemma [5\2|) with respect to t and derive an 
energy estimate for k\. To this end, let us denote 



Then differentiating the curvature equation k\ = k 1 ss + (k*) 3 + v 1 k 1 s with respect to t and 
taking into account the commutation relation dtd s = d s dt we obtain 

w\ = w\ s + 3(k<) V + v l w l s + vy s 

for i = 1,2, 3. Multiplying the above equation with w % and integrating over Y\ yields 



2dt J v 



w l w'ds + 3 I (K i ) 2 (w i ) 2 ds+ I Ww l w\ + v)w l K\)ds 



+ / (v l w l wl + vlw l K l s ) ds . 



wlYds + 3 I (K l Y(w l Yds 



(6.14) 



In what follows, we analyze the boundary term [u^w*] ^ appearing in the right hand 
side of (I6.14p . First we analyze the boundary term at the triple junction position s = 0. 
Differentiating (15.121) with respect to t, we obtain 



5> 

i=i 



w l = 



(6.15) 



at the triple junction point pit). It follows from (15.131) that there exists a function G{t) 
such that 

4(0, t) + k*(0, t)v\0, t) = G(t) for t > and i = 1, 2, 3 . 
Differentiating this equation with respect to t, we conclude 

w\ + wV + k'vI = G'(t) . 

Therefore we obtain, by using (16.151) . 



i=l 



i=l i=l 

3 

-^7 i {(^)V + KW^}. 



(6.16) 



i=l 



By Lemma 16.31 we can express the term v % as a time independent linear combination of 
curvatures k % (i = 1, 2, 3) evaluated at the triple junction pit) and so v\ can be expressed 
as a time independent linear combination of k\ = w % (i = 1,2,3). Therefore there exists 
a constant C > such that 



i=l 



< C||iu||^oo ||k|| L o 



(6.17) 



Next we proceed with the estimation of the boundary term at the point X 1 e T l n 
i.e. we consider s = r l (t). Notice that r % is no longer constant and its dependence on time 
t has to be taken into account. We will differentiate the boundary condition (15.141) 



Ki(r\t),t) + h\X\r\t),t))K\r l (t),t)=0 



with respect to t. Since ^X\r\t),t) = X\ + X^r 1 )' = k l N 1 + {V + (r i )'}T i = tfN 1 (see 
Lemma 15731) and K l sa = k\ — (V) 3 — t>V s =ty ! - (V) 3 — we obtain 

w^ + tiw* = - (k»„ + /iV s ) (r*)' - (VV, 

= {w i - («*) 3 - (h* - v*)/iVy - (j) 2 (Vti, N% 2 . 



(6.18) 



Here we have used the equations + h l n l = and t> * + (r 1 )' = and expressed hi by the 
right hand side of (16.111) . We now denote by S 2 the right hand side of fl6.18p . Then we 

get 

3 3 3 

8=1 i=l i=l 

As the outer boundary dQ is assumed to be C 3 smooth we obtain that the terms h % , Vh l 
are uniformly bounded. Hence the remainder term Yli=i 7 l «/S l can be estimated as 



t=i 



< c\\w\\ L °° (||k|||oo + HHU^IMU 00 + 



+ v 



\f^f\\xjOo -j- ||f ||x°° || ^|| L°° ) ■ 



(6.20) 



In order to complete our estimates we have to derive L°°- estimates on the tangential 
velocity v % and its time derivative v\. Since v\ = {k 1 ) 2 we have 



v\s,t)=v*(0,t) + / |«*(C,*)rdC- 

JO 



(6.21) 



By Lemma 16.31 we can express v l (0,t) as a time independent linear combination of k\ 
(i = 1, 2, 3) evaluated at the triple junction p(t). Therefore there exists a constant C > 
such that 

IMU« < C(\\k\\ L oo + ||/s|||a) < C||«|| L oc (6.22) 

for ||k||l 2 < 1- Analogously, as vl(s,t) = vl(0,t) +2 J s K l ((,t)K l t ((,t) d( and fj(0,t) is 
a time independent linear combination of k\ = w l (i = 1,2,3) evaluated at the triple 
junction position pit), we conclude 



l|i>tlU°° - cdl^lU 00 + II k IU 2 IIHU 2 ) - 

for ||k||z,2 < 1. 

Summarizing we have shown the following equality 



(6.23) 



1 d 
2~dt 



3 

/ (w i ) 2 ds + I[w,w] 

i=i J n 



-«>X +3 / iK i ) 2 (w t ) 2 ds+ / (vWwi + viw^i) ds 



and, consequently, the estimate 

--H| i2 + /[™,™] 

<C{(|| ^|||oo + II^Hlcx; + ||K.||^oo) ||ioo + ||#«||l°o ||lu|||oo + 1 1 K, \ \ \oo || W || £2 

+ll K IU°°ll' u; IU 2 ll' u; slU 2 + II^IU 00 II' u; IIl 2 II k; sIIl 2 } • (6.24) 

The application of the above inequality will be twofold. At first, we utilize it in order to 
prove a bound on ||/t ss (£)||.£,2 uniformly for t G [0, T) where T > is the maximal time of 
existence of a C 2+a solution p. This implies together with Theorem 12.31 the possibility of 
global continuation of the C 2+a solution p up to the maximal time of existence T = +oo 
and hence the global existence of a C 2+a solution will follow. As a second application of 
the above inequality we will prove exponential stabilization of a solution in the if 2 -norm 
of the curvature yielding the exponential stabilization p(t) in its phase-space C 2+Q -norm. 

To accomplish this goal, we have to establish bounds for H^Hl 00 m terms of the norms 
1 1 w 1 1 1,2 and 1 1 k 1 1 £2. This can be done by taking into account the equation w t = n\ = 
K 8s + ( K *) 3 + v l K l s . From this equation we have, for i = 1,2, 3, 

KJl 2 < C(||^|| l2 + ||(^) 3 ||l 2 + ||^<||l 2 ) 

< C(|K|| i2 + ||« i ||£e + IMU-lklU 2 )- (6-25) 

Let us denote by || • ||#fc the following Sobolev norm of the Sobolev space H k = W k ' 2 

HvlU* = II^IU 2 + II^VIU 2 - 

Due to the continuity of embeddings H 2 H 1 and H 2 '—>■ L°° and using Gagliardo- 
Nirenberg interpolation inequalities (cf. [TJ Lemma 5.18 and Theorem 4.17]), we infer the 
existence of a constant Co > such that 

13 11 

|k|| L oo < Co||k.||^- 2 ||k.||| 2 , II K s || L 2 < Co \\k\\^j2 \\ k \\^2, (6 26) 

\w\\ L oo < CoWwW^WwWls . 

By the Young inequality ab < a p /p + b q /q with p = 4/3, q = 4, we have, for any e > 0, 
||«IU°°||«Ji,2 < C 2 ||k||^ 2 ||k||^ 2 < £||k||h2 + C £ ||k||| 2 

and, analogously, ||/«||f,6 < < CoCH^H ^ 2 \\ k \\l2 < ^II^IIh 2 + C e 1 1 1 1 • By taking 

< e <C 1 small enough we obtain from (16.251) 

||« ss ||l2 < C (||k||l 2 + HHU 2 ) for ||k||l2<i. 

Consequently, 

II^IIl 2 < C ||k||| 2 ||«||12 < ^ (IlKllfl, + ||k|U 2 ) < C{\\k\\ l i + \\wWv) 
for ||k||l2 < 1. Similarly 

< C{\\k\\ l i + \\w\\ L 2) 



for 1 1 /« | |z,2 < 1 where C > is a generic positive constant. Due to the continuity of 
embedding H 1 L°° we have ||iu||l°° < C||iu||tfi. 

We proceed by estimating the right hand side of (16.24ft . From ( I6.26f) we have 

||k||1=o <C ||k|| h2 < C(\\k\\ L 2 + \\w\\ L 2) <C{l+\\w\\ L 2) 

for ||k||l 2 < 1- Consequently, by using the Young inequality, we obtain 

||«|||oo + ||«|||oo + ||k||1°o < C (1 + || if || l 2 ) 

for || k\ I £2 < 1. From the embedding ||it>||^oo < C||iu||#i and Young's inequality it follows 
that 

(||«|||oc + ||k||1=c + ||k||jU)||«>||l«> < ^II^IIh 1 + C £ {i + IIHIi 2 ) 

for || At || ^2 < 1. Using the Gagliardo-Nirenberg inequality ( 16. 26ft and Young's inequality, 
we can estimate the second summand in ( 16.241) as 

||k||l°° ll^lli 00 ^ C|MU°° II'HI/HI'^II.l 2 < £ II w \\ m + Cell K lll°° II w \\ z,2. 

Then, by means of ||/c||£oo < H K lli 2 + II^sIIl 2 — 1 + ll^slli 2 f° r ll^-IU 2 — 1> we have 

||K||L«l|w|lia < (! + Kll! 2 )ll™ll! 2 > 

and so 

||w|||,oo < + C £ (i + H^IIl 2 )!!^!!! 2 - 

The remaining terms in (16.241) can be easily estimated with help of Young's inequality as 

IMU^IMUHI^IU 2 < ell^llfl-i + C £ (\ + ||k s ||| 2 )||ii7|||2, 
IIHU^IMUHI^IU 2 < elMlffi + Cell^lll 2 !!^!!!, 2 

for ||k||l 2 < 1. Let us introduce rj(t) := l + ||K s (t)||| 2 - Then, by choosing e > sufficiently 
small and taking into account the positivity of the bilinear form /, we obtain 

I[w,w] > S\\w\\ 2 H i. 

Therefore the function ||iy(t)||| 2 satisfies the differential inequality 

^IM*)H£» <C 1 + C 2V (t)\\wml 2 . (6.27) 
According to Lemma MM the function rj is integrable on the interval (0,T) and 

pT p+oo 

/ Tf(t)dt<T+ \\K 8 (t)\\ 2 L 2 dt < oo 
Jo Jo 

provided that T < oo. A Gronwall lemma type of argument applied to the differential 
inequality (16.271) yields the existence of a CV, which is monotone increasing and bounded 
as long as T is bounded, such that 

sup \\w(t)\\ 2 L 2 < C T < +oc . 

0<t<T 

By means of Lemma [6.61 we see that ||k||l2 is small if ||p(0) ||c* 2 + a is small enough. Fur- 
thermore, ||p||c 1 + a (0 < ot < 1/2) is small provided when ||k||l 2 is small. In addition, 
using ||k ss ||^2 < C (||k||l 2 + ||ty||L 2 ) and the fact that the norm ||p||co+a can be estimated 
by ||k||jj2, we just have shown the following conclusion. 

Theorem 6.7 The local solution of Theorem \2.3\ can be extended to the time interval 
[0, oo) provided that p is small enough in the C 2+a -norm. 



7 Exponential stability of stationary solutions 

In this section we combine all the previous results to prove exponential stabilization of a 
solution to the triple junction problem which have initial data close to a stationary stable 
solution. 



Theorem 7.1 Let the assumptions of Theorem \2.0A and Theorem\6.1\ hold and let be 



such that the bilinear form 7* is positive. Then there exist constants C , uo , S > such 
that 

\\p(t)\\ H 2 <Ce- wt ||K(0)|U 2 
for any t > and \\k(0)\\l3 < 5. 

Proof. The proof directly follows from Lemma 16.61 and Corollary 14.21 □ 

Since the i? 2 -norm of p dominates its C 1+Q -norm and the C 2+Q -norm majorizes L 2 - 
norm of k we can state the following consequence of the previous theorem. 



Corollary 7.2 Under the assumptions of Theorem \7.1\ there exist constants C, u,5 > 
such that 

||p(t)|| C i +a <Ce-"'||p(0)|| C2+ « 
for any t > and ||p(0) ||c ,2 +« < 8- 

Finally, we are able to prove exponential decay in stronger norms. As it was already 
indicated in the previous section, we will utilize the higher energy estimate (16.241) once 
more in order to prove exponential stabilization in the if 2 -norm of the curvature k. 

Recall that, for p > 2, we have 

H K IISollHU°° ^ c p (IMU 2 + IIHU 2 ) P IIHU°° 

< c p (|H| L 2 + WwW^p-WIkW^WwWl^ + |H| L2 |H| L oo) (7.1) 

< C(||k|| l2 + ||iu|| L 2)(||K||| 2 ||if;||^ 1 ) 

< c\\k\\ 2 L 2 + c(||k|| L 2 + IMMIMItfi 

provided that || k ||^2 < 1 and ||iw||l2 < 1. Since || w \\ L 2\\w s \\ ip. < H^H^i and || || 1 1 it> 1 1 ^2 < 
CIHI 2 ^ we conclude from (1fT24l) . <K2~Et\ . (1731) 

~\\w\\ 2 L2 + I[w,w] < C\\k\\ 2 l2 + C(\\k\\ L 2 + ||w||La)||tw||^ 

for some positive constant provided that ||k||l2 < 1 and ||iu||l2 < 1. 

Similarly as in the proof of exponential decay of ||k||l 2 we use the fact that the full 
Sobolev norm ||iu||#i can be estimated by the bilinear form I(w,w) as follows: 

^II^II/p < I[w, w] 

for some positive constant 5 > 0. Taking ||k||l 2 and ||iu||l 2 sufficiently small such that 
C ( 1 1 ^ 1 1 ^,2 + || it7 1| z,2 ) < 5/2 we end up with the inequality 



1 d n 112 5 .. , 1 



Defining y = \\w\\^ 2 and using ||iu||l 2 < ll^llir 1 we have 

rl.ii 

+ Sy< 2CMe- ut 



dt 

where M, u > are the modulus and rate of exponential decay of ||k||^ 2 , i-e. || /«.(., t)\\ 2 L 2 < 
Me" ul . Solving the above differential inequality with respect to y = y(t) we end up with 
the following estimate 

y(t)<y(0)e-* + ^L\e-^-e- St \ . 

It means that the norm ||io||^2 exponentially decays with the rate min(5, uj). Since 
|| k ss||l 2 < C(II K IU 2 + II^IU 2 ) and the full Sobolev norm H^Hh 1 dominates as 
well as the ||K s || L 2-norm we obtain the following convergence result: 



Theorem 7.3 There exist constants S, M, uj > such that the solution of Theorem \6.1 
fulfills 

\\H,(t)\\ 2 H2 < Me- Wt 
for allt>0 provided that ||k(0)||^- 2 < S. 
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